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Abstract 

We study elementary eigenfunctions y = pe^^ of operators L{y) = 
y" + Py, where p, h and P are polynomials in one variable. For the 
case when h is an odd cubic polynomial, we investigate the real level 
crossing points and asymptotics of eigenvalues. This study leads to 
an interesting identity with elementary integrals. 
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1 Introduction 

Following Bender and Boettcher [3], we consider the eigenvalue problem 

iv" + {C^ + 2bC + 2tJC + X)w = 0, w;(te-*("/2±"/^)) ^ 0, t ^ +oo, (1) 
*Both authors were supported by NSF grant DMS- 1067886. 
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with real parameters b and J. 

When J is a positive integer, this problem is quasi-exactly solvable (QES) 
which means that for every b there are J eigenvalues Xj^k, k = 1, . . . , J with 
elementary eigenfunctions 

yj,kiC) = PjAC) exp(-<V3 - ibC), 

where pj^k are polynomials of degree J — 1. These J eigenvalues are found 
from algebraic equations 

Qj(6,A) = 0, J=l,2,..., (2) 

where Qj are real polynomials of degree J in A. 

For every J and b, problem ([T]) has infinitely many eigenvalues A tending 
to infinity, which satisfy an equation 

Fj{b,X)=0, (3) 

where Fj is a real entire function of two variables which is called the spectral 
determinant. 

We denote by Zj{lV) and Zj^'^(R) the sets of real solutions (&, A) of 
equations ([3]) and ([2]), respectively. These are certain analytic curves in R^. 
A computer-generated picture of these curves is given in [3]. The symbols 
Zj and Z^^^ denote the sets of all complex solutions of these equations. 

The most conspicuous feature seen in the picture in ^ is a set of real 
crossing points where the curves of Zj^'^(R) intersect the rest of Zj(R). 

In this paper, we prove that infinitely many such points exist for each 
odd J (Proposition 1). We characterize real and complex crossing points as 
intersection points of z'j^^ with Z_j (Theorem 1). We prove that all non- 
QES eigenvalues are real when b is real and J is a positive integer (Theorem 
2.) 

In sections [H [9], we study the asymptotic behavior of polynomials Qj as 
b — )■ ±oo, using singular perturbation methods from [TO]. Appendix A with 
a description of the relevant results of [10] is included. 

Our proof of Theorem 1 is related to a challenging conjecture on certain 
elementary integrals. We rigorously verified this conjecture for integers J < 5 
using symbolic computation. The conjecture is discussed in sections [3l [51 
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2 Elementary solutions 



Let h and p be two polynomials in one variable. When y = p{z)e^^^'^ satisfies 
a second order differential equation 

y" + Py = \y, (4) 

where P is a polynomial? Substitution gives 

^ + 2^h' + h" + h'^ + P- X = 0. (5) 
p p 

Such P exists if and only if 

p" + 2ph' is divisible by p. (6) 

Another criterion is obtained if we consider the second solution yi of (j4]) 
which is linearly independent of y. This second solution can be found from 
the condition 

yy'i - y'yi = i- (7) 

Solving ([7]) with respect to yi we obtain 

y,=pe^ jp-'e-^\ (8) 

As all solutions of ^ must be entire functions, we conclude that 

all residues of p~'^e^^^ vanish. (9) 

This condition is necessary and sufficient for y = pe^ to satisfy equation 
dl]) with some P. Indeed, if (j9]) holds, then yi defined by (j8]) is an entire 
function, so {y, yi) is a pair of entire functions whose Wronski determinant 
equals 1, so this pair must satisfy a differential equation (jlj) with entire P, 
and asymptotics at infinity show that P must be a polynomial. 

Thus conditions and ([9]) are equivalent. One can give another equiv- 
alent condition in terms of zeros of p, as in [19] . Let 

n 

P{z) = Y[{z-Zj), Pk{z) = p{z)/{z - Zk). 

J=l 

Then (jG]) is equivalent to 

p"{zk) + 2p'{zk)h'{zk) = 0, 
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for all = 1, . . . , n. We have p'{zk) = Pk{zk) and p"{zk) = 2p'^(zfc), so the 
condition 

E = -h'{zk), l<k<n, (10) 

is equivalent to ([6]) and ([9]). Equation (|TOl) is the equilibrium condition for n 
unit charges at the points Zk in the plane, repelling each other with the force 
inversely proportional to the distance, and in the presence of external field 
h'{z). Equations ffTOj) express the fact that {zi, . . . , Zn) is a critical point of 
the "master function" 

^{z,,...,Zn)= n {z,-z,)Jle^^^-\ 

{j,k):k<j k 

3 Remarkable identity 

From now on we suppose that h is an odd polynomial of degree 3, which we 
write in the form 

h{z) = zy3-bz. (11) 

Suppose that all residues of p~'^e^'^^ vanish. Then the integral / p~'^e~'^^ is a 
meromorphic function in the plane. Surprisingly, the integral of some linear 
combination 

j (p2(_;,)e-2M^)_Cp-2(^)e-2M-)) 

is not only meromorphic but is an elementary functionl Here C is a constant 
depending on b and p. 

Conjecture. Let h be given by Let p be a polynomial. All residues 

of p^^e"^'* vanish if and only if there exist a constant C and a polynomial q 
such that 

In other words: 

p'^{z)p'^{—z) — C = q\z)p{z) — q{z)p'{z) — 2q{z)p{z)h'{z). 

It is known [3] that for given h of the form ffTTl) there exist polynomials p of 
any given degree such that all residues of p~'^e~'^^ vanish. These polynomials 
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p have simple roots. We verified the conjecture for degp < 4 by symbolic 
computation with Maple. We don't know whether there is any analog of the 
Conjecture for other polynomials h. 

Substituting Pn{z) = z"' + az"'~^ + . . . into ([5]) and using ffTTl) . we conclude 
that 

P{z)-\ = -h'^{z)-h"{z)-2nz + 2a = -z* + 2z%-2{n+l)z-h^ + 2a. (13) 

We choose A = 6^ — 2a so that P(0) = 0. Equation (E]) now becomes 

K + 2{z^ - h)p[, - (2nz - 2a)pn = 0. (14) 

Coefficients of pn can be now determined by a linear recurrence formula. 
Putting 

n 

Pn{z) = "^ttjZ^'^, a_i = 0, ttQ = 1, ai = a, 

j=0 

we obtain the recurrence 

h( -^9^ , (n-j + 2)(n-j + 3) 
jaj = aaj_i - b{n - j + 2)aj_2 H flj-s- (15) 

Coefficients aj are found from this formula one by one beginning from ci = 
a. Vanishing of the constant term in (fT4|) gives a polynomial equation 
Qn+ii^j a) = in which we can substitute a = (6^ — A)/2 and write it as 

g„+i(6,A) = 0. (16) 

Polynomials Pn, Qn+i small n are given in Appendix B. 

We have deg^^Qn+i = n + 1, [3]. For every b and every A satisfying this 
equation, the differential equation @, with P as in f|T3|) . has a unique solu- 
tion y = pne^ where pn is a monic polynomial of degree n. Coefficients of pn 
are polynomials in h and A. 

Polynomials Qn+i are fundamental for our subject, but little is known 
about them. It seems hard to investigate them algebraically. In section |H1 
we will use analytic tools to establish some properties of these polynomials, 
in particular we will find the terms of top weight and asymptotics of A as 
6 — )■ oo. 

Functions y = Pn^^ are eigenfunctions of the operator 

Lj{y) = y" - {z^ - 2hz' + 2Jz)y, J = n + 1. (17) 
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This operator maps the space {pe^ : degp < n} of dimension n + 1 into 
itself. For each non-negative integer n, and generic b, the operator (1171) has 
n + 1 eigenf unctions of the form p„e'^ with eigenvalues A which are solutions 
of ([I6D. 

We assume without loss of generality that Qn+i is monic as a polyno- 
mial in A, and p„ is a monic polynomial in z. Then the constant C in the 
Conjecture turns out to be 

C{b,X) = a^^Q^+,. (18) 
Symbolic computation for small n shows that a„ = (— 1)"2^^". 



4 Boundary value problem 

Eigenf unctions pe^ do not belong to L^(R), but they satisfy the boundary 
conditions 

yite^^''^) -^0 as t -> oo. (19) 

With these boundary conditions, the operator ( IT71) is not Hermitian but PT- 
symmetric |3l HH |18] . 

Physicists write the boundary value problem for the operator Lj with 
the boundary conditions f lT9|) in the equivalent form ([1]), which corresponds 
to the rotation of the independent variable iC, = z, w{() = y{i()- We also 
find this form convenient in certain arguments, and will use it in sections [TJ 
m We keep the notation y{z) for an eigenfunction of ( !T7|) . (fT9|) . while w{() 
stands for an eigenfunction of ([1]). 

It is known that the boundary value problem ( IT7|) . ( IT9|) has an infinite 
sequence of eigenvalues tending to infinity [16]. Eigenvalues A are solutions 
of the equation 

F„+i(6,A) = 0, (20) 

where is a real entire function on which is called the spectral deter- 
minant [13]. The set of all solutions of (120|) in is called the spectral locus 
and we denote it by Zn+i- As Fn+i is real, the set of eigenvalues is symmetric 
with respect to the real line when b is real. For each real b, all sufficiently 
large eigenvalues (how large, depends on b) are real |14] . 

Equation ( 120|) is reducible: is evidently divisible by Qn+i- On the 

other hand, equation f|T6|) is irreducible, as follows from [9] or [2], and it 
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defines a smooth algebraic curve in C^. This algebraic curve will be denoted 

by z^^^t 

5 Theorem 1 

Now we discuss a corollary of our conjecture that we can prove. Let us fix 
a simple curve 7 in C parametrized by the real line, with the properties 
7(t) -7- 00, arg7(t) — > ±7r/3, t — > ±00. Then f lT2l) implies 

p\z)e^^^'^dz = C [ p-\-z)e^^^'Uz. (21) 
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To obtain this we replace z 1— )■ —z in f[T^ then integrate along 7; the integral 
in the right hand side of f[T^ vanishes because ^h{z) — —00 as z — )■ 00 on 
7. Let 7^ be a curve consisting of the piece {'jit) : —00 < t < 0} followed by 
a curve from 7(0) to z. Put 

(^(z) = p(_2)e-^(-) / p-2(-C)e2'^(^)rfC- 
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Then g^{z) := g{—z) satisfies Ln+i{g^) = Xg" so g satisfies L_n-i{g) = 
Xg. To check this, we make a substitution z ^ —z in ( 1T7|) . Moreover, 
if the integral in the right hand side of ( 1211) is zero, then g also satisfies 
the boundary condition ( !T9|l . Therefore our conjecture ( !T2|) . ( !T8|) has the 
following corollary: 

Theorem 1. The points {b, X) G Z^^i where the eigenfunction y = pe^ 
satisfies 

''y^{z)dz = (22) 

are either zeros ofC{b,X) or points of intersection of Z^+f with Z_n-i- 
We will prove Theorem 1 in the next section. 

Equation fl22|) is the well-known condition of level crossing, which we 
discuss in section [71 

Thus the Corollary says that the eigenvalues at the points on Z^^i which 
are singular points of ^„+i are eigenvalues of two spectral problems, one for 
Ln+i another for 
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6 Darboux transform 

Proof of Theorem 1. Assuming that all residues of p~^e~^^ vanish, we will 
prove that the right and left sides of (!2T|) with C as in (fT8|) have the same 
zeros on Z^j^i . Fix the integer n > 0. Let ipk = PfcC , k = 0, ... ,n be 
all elementary eigenfunctions of Ln+i- They are linearly independent, and 
they span a space V invariant under As \^ is a subspace of U = 

{pe^ : degp < n}, we conclude that V = U. So the Wronski determinant 
W = W{ipo, . . . , ipn) is proportional to the Wronski determinant 



Now let us perform the Darboux transform of Ln+i killing these n+1 eigen- 
functions. We recall that Darboux transform [5l [131 IH E]) applies to any 
operator —D"^ + V with eigenfunctions 0o, • • • , 0n and corresponding eigen- 
values Aq, . . . , A„. The transformed operator is 



and its eigenvalues are those eigenvalues of —D^ + V which are distinct 
from Ao, . . . , A„. As in our case 2{\ogW)" = 2{n + \)h" = 4:{n + l)z, the 
result of application of the Darboux transform to L^^i and eigenfunctions 
ipk, k = 0, . . . ,n is L_„_i. 

If the left hand side of ( l2Ti) is zero at some point (6, A) G Z^+f then 
dFn+i/d\ = at this point (see, for example 03]). As Z^+i is smooth, this 
is possible in exactly two cases: either (6, A) is a smooth point of Zn+i and 
dQn+i{,b, X)dX = 0, or (6, A) is a self-intersection point of Zn+i- 

In the second case, {b, A) belongs to the spectral locus of the Darboux 
transform L_„_i. This means that the equation 



with these parameters (6, A) has a solution y* that tends to on both ends of 
7. Then yi = y*{—z) tends to on both ends of —7 and satisfies L„+i(?/i) = 
Xyi. So yi satisfies the same differential equation Ln+i{y) = Ay as y does, 
and is linearly independent of y. So yi = y J y~'^e~^^. As this tends to on 
both ends of —7, we conclude that / y~'^e^^^ tends to on both ends of —7. 




Z}2 + y_2— logiy(0o,...,0n), 



L 



— n— 1 



(y*) = Ay* 
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So y*{z) = yi{—z) tends to on both ends of 7, and this means that the 
right hand side of (121 p is 0. 

The argument is evidently reversible. This proves that the right and 
the left hand sides of f l2T|) have the same zeros on Z^^^ , that is f l2T|) with 
C = andQn/dX, where A) 7^ on Z^+i- 

According to the theorem of Shin pTS], all eigenvalues of Lj for J < 
are real. Shin's proof of this uses the ODE-IM correspondence discovered by 
Dorey, Dunning and Tateo, [B]. 

Combining the Darboux transform used in the prof of Theorem 1 with 
the result of Shin [15], we obtain 

Theorem 2. For every positive integer J, all non-QES eigenvalues of Lj 
with boundary conditions / fT^) are real. 

Proof. These eigenvalues are also eigenvalues of L_j with boundary con- 
ditions f liyp . and the eigenvalues of L_n-i are all real by Shin's theorem [T3] . 

For J > 1 there are always some non-real eigenvalues. 

7 Level crossing 

As Qn+i and Fn+i are real functions, it is reasonable to consider real solutions 
of equations ( 120|) and (fT6|) . Eigenfunctions y{z) corresponding to these real 
solutions are real, while eigenfunctions w{() (see ([1])) are PT-symmetric, that 
is w{—() = w{(). These real solutions (6, A) form curves in which we call 
the real spectral locus Zn+i(R) and the QES real spectral locus Zn+i(R), 
respectively. 

Now we discuss fl22l) . First we state a result which describes Z^^f(R). 

Theorem 3. For n > 0, the spectral locus Z^^f(R) consists of [n/2] + 1 
disjoint analytic curves Tnm, < m < [n/2] (analytic emheddings of Ti to 

For {b, A) G r„ the eigenfunction has n zeros, n — 2m of them real. 

If n is odd then b — )■ +00 on both ends of each curve Tn,m- If ^ is even 
then the same holds for < m < n/2, but on the ends of Tn,n/2 we have 
b ±00. 

// (6, A) G Tn,m^ {b, ^) G r„^m+i o,nd b is sufficiently large, then > A. 
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The proof of this theorem can be found in [TT]. It follows the method 
of [TU] where similar results were established for real spectral loci of other 
families of cubic and quartic potentials. The method is based on singular 
perturbation and Nevanlinna parametrization of the spectral locus. 

Computer generated pictures of Z„+i(R) show an interesting phenome- 
non: when n is even, the curve r„ „/2 crosses the non-QES part of the spectral 
locus Pl Fig. 1]. We will prove that infinitely many such crossings exist for 
even n and negative h. 

We say that a level crossing occurs at a point (6, A) of the spectral locus 
if dFn+i/dX = at this point. If y is the eigenfunction corresponding to a 
point (6, A), then the level crossing occurs if and only if (12^ is satisfied 
II. 7], fIE[ Thm. 8]. There are two types of level crossing points: 

a) Critical points of the function A at non-singular points of Zn+i- 

If such a critical point {bo, Aq) is simple and belongs to Z„+i(R) then the 
two eigenvalues that meet at this point are both real for b on one side of bo 
and complex conjugate on the other side. 

b) Singular points of Z^+i- 

If two eigenvalues collide at a simple self-intersection point of Z„_|_i(R) 
with two distinct non-vertical tangents, then these eigenvalues both remain 
real in a neighborhood of bo. Operator Ln+i with b = bo contains a Jordan 
cell in this case. 

We recall that Z^^f is a smooth curve. Thus the crossing points on Z^^f 
where only QES eigenvalues collide are all of type a), and they satisfy 



For each n, there are finitely many such points on Z^+i . 

We will show that there are always infinitely many crossing points of type 
b) where QES eigenvalues collide with non-QES eigenvalues. So the curve 
defined by fl2Up is not smooth: it has infinitely many self-intersections. 

We don't know whether more complicated singularities than a) and b) 
exist; numerical experiments only show singularities of types a) and b). 

Proposition 1. Function 



Q„+i(6,A) = 



dX 



Q„+i(6,A) =0. 
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where y is the eigenfunction corresponding to {b, A), has infinitely many zeros 
{bk, Afc), bk — > oo. When n is even, has infinitely many zeros with negative 
bk and real Xk- 

Proof. We have 

$n(&) = / pl{z)e'''^^^dz. 

We remind that coefficients of p„ and h are algebraic functions of b. When 
n = 0, we can take po = 1, and then 

$0(6) = / e^^^^^^'-^'^dz = 22/3^7rAi(22/36), 

where Ai is the Airy function pj. Airy function is a real entire function of 
order 3/2 with infinitely many negative simple zeros. 

To generalize this to other values of n, we express as a linear combi- 
nation of $0 and $0 "with coefficients depending on b algebraically. Differen- 
tiating $o(&) with respect to b, we obtain 

and thus 

$„(6)=p^(-D/2)$o(6), 

where D = d/db. Now all ^^q^ are linear combinations of $0 and $q with 
polynomial coefficients because Ai satisfies the differential equation Ai"(s) = 
sAi(s). So is of the form 

<1>„(6) = A„(6)<l>o(6) + Br,m',{b), (23) 

where An and Bn are algebraic functions. 

We claim that every linear combination of $0 and with algebraic 
coefficients has infinitely many zeros. We prove this claim by contradiction. 
Suppose that such a linear combination 

(f) = ao$o + ai'^'o (24) 

has finitely many zeros. Let F be a compact Riemann surface spread on the 
Riemann sphere on which oq and ai are meromorphic. Then (p is meromorphic 
on F\E, where E is the finite set of points of F lying over 00. At the points 
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of E, (f) has isolated essential singularities. As has finitely many zeros 
and poles on F\E, we conclude that (f)'/4> is meromorphic on F\E. The 
growth estimate log|0(6)| < 0(|6|'^/^), b — )■ oo, implies that the points of E 
are removable singularities of (^'/(p. Thus is the exponent of an Abelian 
integral. Now consider fl2^ as a linear differential equation of first order with 
respect to $0) whose coefficients belong to the minimal field K that contains 
C(6), is algebraically closed, and contains a primitive of every element, and 
the exponent of a primitive of every element. As every first order linear 
differential equation can be solved by integration we conclude that E K 
which implies that Ai G K. But this is not so by a well-known classical 
theorem of Picard and Vessiot, [121 Theorem 6.6]. This proves our claim. 

When n is even, according to Theorem 3, we have a real analytic branch 
A (6) defined for all real b with sufficiently large absolute value. The graph 
of this branch is a part of r„^„/2- Using this branch we rewrite the equation 
$n(6) = as 

$'o(6)/$o(6)=A(6), 

where A is a real branch of an algebraic function on {—oo,B) with some 
-B G R. This last equation has infinitely many negative solutions because $o 
has infinitely many negative zeros and they are interlaced with zeros of $q. 
This completes the proof of the proposition. 

Using the asymptotics of the zeros of Airy's function [1] we obtain that 
the crossing points satisfy b^ ~ — ((3/4)7rA;)^/^, A; — > oo. 

8 Asymptotics as 6 — > +00 

Now we study asymptotics of the eigenvalues A as 6 — ?■ +00 and make conclu- 
sions about polynomials Qn+i- Our main result here is the explicit formula 
fl5U]l for the top quasi- homogeneous part of Qn^i- 

First we obtain a preliminary estimate of solutions A (6) of equation ( !T6|) 
for large b: 

X{b) r^b^ + 0{Vb), 6->oo. (25) 

To prove this, consider the recurrence ffTSj) . For a monomial a"^b^ we define 
the weight as m -|- 2k. Then ( lT5l) implies that 

[i/2] 

jlaj = + ^ Cm,jb'^a^~'^"^ + terms of lower weight. 

m=l 
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Vanishing of the constant term in fll4p gives 

(5*+i(a, b) = aan + &a„_i + 2a„_2 = 0, 

so Qn+i is a sum of a quasi- homogeneous polynomial in a and b of weight 
2(n + 1) and a polynomial of lower weight. This means that a = 0{\/b) and 
A (6) = b^-2a satisfies 025]). 

To obtain more precise asymptotics we use singular perturbation argu- 
ments from [lOj, which we state in Appendix A for the reader's convenience. 

Suppose that b is real and b — i- +oo. In the equation ([1]) we set 

C = CM — ze^^, 6 = e~^, W{u) = w{eu — ie^^). 

The result is 

W" + (eV - Aie^u^ - ^v? - 2iJe\)W + (2J + e^A - e-^)W = 0, (26) 

or 

- W" - (^^(r^/^M - 2^)2 - 2iJb~^/^u) W = {2J + r^/'A - b''/'')W. (27) 

When e — 7- 0, we obtain the limit eigenvalue problem 

- W" + Au^W = nW, (28) 

which is a harmonic oscillator with eigenvalues = 2(2A;+1), k = 0,1,2, .. .. 
By a general result from [lOj (see Appendix), fl27|) implies that for each k, 
there must be a unique eigenvalue Xk{b) which satisfies 

Xk = b^ + {l2k - 2J + o{l))Vb. (29) 

Moreover, for each compact set K in the A-plane there exists b^ > such that 
for b > bo there are no other eigenvalues A (6) G K, except those satisfying 

We conclude from ( |25ll that QES eigenvalues must satisfy ( l29l) . That is 
for each QES eigenvalue A there exists k such that fl29l) holds. Now we have 
to find out what are the values of k for the QES eigenvalues. 

To do this, we consider zeros of eigenf unctions. We know that k-th eigen- 
function of fl28l) has [k/2] zeros in the right half-plane, the same number of 
zeros in the left half-plane, and one zero on iR if k is odd. (In fact all these 
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last zeros belong to the real line but this is irrelevant for our argument.) So 
for every m = 0, 1, . . . there are two eigenf unctions of the harmonic oscillator 
(with k = 2m and k = 2m + 1) which have m zeros in the right half-plane, 
and one of them {k = 2m + 1) has a zero on zR. 

Theorem 3 implies that for each given n and for each m < [n/2] and 
b sufficiently large positive, there is exactly one curve Tn^m, such that the 
corresponding eigenfunctions have m zeros in the right half-plan^. We refer 
to [TOj for the argument showing that the zeros of eigenfunctions w in the 
right half-plane do not escape to infinity as 6 — )• +oo. Zeros of w on iH do 
escape to infinity, except possibly one of them. Thus the branches of QES 
eigenvalues must be Aq, . . . , A„ satisfying (12^ . 

Putting Afc = 6^ — 2a{k), and J = n + 1 in f l29p we obtain 

a{k) Vb{n-2k), < k < n. 
We conclude that the top weight term of the polynomial Q^+i is 



f\ (a-(n- 2k)Vb] =1 ~ ^^^^^ - 36) . . . (a^ - n6), n is odd, 
V ^ I — 26) . . . (a^ — n6), n is even. 

(30) 

This implies that the degree of the discriminant of Q^+i is n(n+ 1)/2, and the 
genus of the QES spectral locus is n{n — 2)/4 when n is even and {n — 1)^/4 
when n is odd. 



9 Asymptotics as 6 — > — oo 

When b — > — oo, our operator ([T]) also degenerates to a harmonic oscillator. 
However none of the QES eigenvalues of ([T]) tend to the eigenvalues of this 
harmonic oscillator as 6 — )■ — oo. To study this limit, we set z = eu, b = — 
and W{u) = w{eu) in ([1]). The result is 

W" + [e\^ - 2u^ + 2iJe'u + e^A) W = (31) 
As e — > 0, this tends to the harmonic oscillator 

-W" + 2u^W = nW, 

^Remember that we are working here with eigenfunctions w{C,) ~ y{iC), where y is an 
eigenfunction from Theorem 3. 
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whose eigenvalues are fXk = \/2{2k + 1), A; = 0, 1, 2 ... . So by the results in 
[To] (see Appendix A), for every k and for b < —bk, there is an eigenvalue 
Xk{b) which satisfies 

\ime'^Xk{b) = V2{2k+1), 

or Xk{b) ~ ^/^. Comparison with fl25|) shows that these eigenvalues 
cannot come from the QES spectrum. 

We thank Per Alexandersson for making Fig. 1, and for help with com- 
putations which led us to the discovery of flT^ . f lTS]) . Stefan Boettcher for 
sending us pictures of Zj(R) which inspired this work, Evgenii Mukhin and 
Alexandre Varchenko for useful discussions, and Vladimir Marchenko for his 
insightful suggestion to look at the Darboux transform. 
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Appendix A. Singular perturbation of poly- 
nomial potentials 



Here we state the main singular perturbation result of |TU] and verify 
that the eigenvalue problems fl26l) and fl3T]) satisfy all conditions that imply 
continuity of the discrete spectrum at e = 0. 

Consider the eigenvalue problem 

-y" + P,{z,h)y = \y, y{z) ^ zeRiUR^. (32) 

Here z is the independent variable, P is a polynomial in z whose coefficients 
depend on parameters e > and 6 G C, dependence on b is holomorphic, 
and Ri, R2 are two rays in the complex plane defined by Rk = {te^^'^ G : 
t > 0}, = 1,2. 
Suppose that 

d 

j=0 

where arf(e) > does not depend on b, Po{z,b) = am{b,0)z"^ + . . ., where 
m < d, and the dots stand for the terms of smaller degree in z. 
Let ^ 

P:{z,b)=Y.ci,{b,e)z^. 

j=m 

For every polynomial potential P{z) = anZ^ + ... of degree n, the sepa- 
ration rays are defined by 

{zeC: anz""^^ < 0}. 

Turning points are just zeros of the potential P in the complex planeJl Ver- 
tical line at a point z is the line defined by P{z)dz'^ < 0. If P depends on 
parameters, then the separation rays, turning points and the vertical line 
field depend on the same parameters. 

We assume that there exists 6 > and eo > and a compact K C Cb, 
such that for all e G (0,eo) and for all b E K and k G {1,2} the following 
conditions are satisfied: 

(i) I argz — > 5 for all turning points z G C\{0} of P*, 

^This terminology is somewhat unusual but convenient here. In the standard termi- 
nology turning points are zeros of P — A. 
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(ii) For every point z & R^, the smallest angle between and the vertical 
line with respect to P* at this point is at least 6. 

(iii) Rk are not separation rays, for P^, e > or Pq. 

(iv) All coefficients aj{b,e) are bounded from above and \am{b,e)\ is 
bounded from below. 

Theorem A. If the conditions (i)-(iv) are satisfied, then the spectral deter- 
minant of the eigenvalue problem converges as e ^ to the spectral 
determinant of [3^) with e = 0; 

^ Fo, e ^ 0, 
uniformly for {b, X) E K x Ki, for every compact Ki G Cz- 




Fig 1. Stokes complex of P*. 



Now we verify that the family of potentials in (l26l) satisfies all conditions 
(i)-(iv) with d = 4, m = 2. We have 

P*{z) = -e^/ + Aie^z^ + 4z\ 

Pq*(z) = 4z^. The turning points are and 2ie~^. The separation rays are 
aigz G {0, vr, ±7r/3, ±27r/3} for P*, e > (shown in thin sohd lines in Fig. 1), 
and aigz G {±7r/4, ±371/4} for Pq (dashed lines in Fig. 1). The normalization 
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rays are arg^; G {— 7r/2 ± 7r/3}. The bold lines in Fig. 1 represent the Stokes 
complex, that is the integral curves of the vertical direction field P*{z)dz'^ < 
that are adjacent to the turning points. 

Thus conditions (i),(iii) and (iv) evidently hold. It remains to verity (ii). 

To do this we parametrize Ri as z = te~^'^^^ : t > and find the 
direction of the line field aigdz at z by inserting this parametrization to 
a.Tg{P*{z)dz'^) = TT. We obtain 

aig P;{z) e (-7r/2,7r/3), ±aTgdz^ E (27r/3, 47r/3), 

so the angle between dz and Ri is at least tt/Q. Verification for R2 is similar. 
We leave to the reader to verify that conditions of Theorem A are satisfied 

for mi. 



Appendix B. Explicit expressions 

We remind that A = 6^ — 2a where Pn{z) = + az"''^ + . . .. Since 
Qn+i{b, A) and Q*^_^_i{b, a) are normalized so that they are monic polynomials 
in A and a, respectively, we have Q*{b,a) = (-l)"+^Q„+i(6, 6^ - 2a)/2"+^ 
We use the notation C*{b,a) = C{b,b'^ — 2a), where C is the constant from 
f[T^ . Then flTSl) can be rewritten as 

C*(6,a)=<^Q:+i, (33) 

where < = (-l)"2"a„ = 2"". 

Here are results of symbolic computations with Maple. 
For n = 1: 

Pi{z) = z + a, 
Ql{b, a) = a? — b, 
C*{b,a) = a. 

For n = 2: 




Aab + 2, 
4 da ' 
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For n = 



3: 



P2[Z) 



z +az + { -a b] z ab + -a + 1 



7 



1 



6 



6 



Ql{b, a) = - lOa^b + 12a + 96^ 



C*{b,a) 



1 3 5^3 IdQl 
-cb^ ab + - = — 

2 2 2 8 aa 



For n = 4: 



1 



1 



z^-az-^+ { -a-2b] z-{2 + -a-^ ab z--a^b+b^ + -a-\ a 



6 



1 



4 24 



Ql(b, a) = A2a^ - 966 - 20a^b + 64a6^ + 



C*(6,a) 



21 



1 dQt 



—a + —a^ - —a^b + Ab^ — - 

4 16 4 16 da 
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